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ABSTRACT 

Since the classical work by Purcell (1979) it has been generally accepted that 
most interstellar grains rotate suprathermally. Suprathermally rotating grains 
would be nearly perfectly aligned with the magnetic field by paramagnetic 
dissipation if not for "crossovers" , intervals of low angular velocity resulting 
from reversals of the torques responsible for suprathermal rotation; during 
crossovers grains are susceptible to disalignment by random impulses. 

Lazarian and Draine (1997) identified thermal fluctuations within grain 
material as an important component of crossover dynamics. For a £ 10 _5 cm 
grains, these fluctuations ensure good correlation of angular momentum before 
and after crossover resulting in good alignment, in accord with observations of 
starlight polarization. In the present paper we discuss two new processes which 
are important for the dynamics of a £ 10~ 5 cm grains. The first - "thermal 
flipping" - offers a way for small grains to bypass the period of greatly reduced 
angular momentum which would otherwise take place during a crossover, 
thereby enhancing the alignment of small grains. The second effect - "thermal 
trapping" - arises when thermal flipping becomes rapid enough to prevent the 
systematic torques from driving the grain to suprathermal rotation. This effect 
acts to reduce the alignment of small grains. 

The observed variation of grain alignment with grain size would then 
result from a combination of the thermal flipping process - which suppresses 
suprathermal rotation of small grains - and due to H 2 formation and starlight - 
which drive large grains to suprathermal rotation rates. 

Subject headings: ISM: Magnetic field, Dust, Extinction - Polarization 
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1. Introduction 

One of the essential features of grain dynamics in the diffuse interstellar medium 
(ISM) is suprathermal rotation (Purcell 1975, 1979) resulting from systematic torques that 
act on grains. Purcell (1979, henceforth P79) identified three separate systematic torque 
mechanisms: inelastic scattering of impinging atoms when gas and grain temperatures 
differ, photoelectric emission, and H2 formation on grain surfaces^]. The latter was shown 
to dominate the other two for typical conditions in the diffuse ISM (P79). The existence of 
systematic H 2 torques is expected due to the random distribution over the grain surface 
of catalytic sites of H2 formation, since each active site acts as a minute thruster emitting 
newly-formed H2 molecules. 

The arguments of P79 in favor of suprathermal rotation were so clear and compelling 
that other researchers were immediately convinced that interstellar grains in diffuse 
clouds should rotate suprathermally. Purcell's discovery was of immense importance for 
grain alignment. Suprathermally rotating grains remain subject to gradual alignment 
by paramagnetic dissipation (Davis & Greenstein 1951), but due to their large angular 
momentum are essentially immune to disalignment by collisions with gas atoms. 

Spitzer & McGlynn (1979, henceforth SM79) showed that suprathermally rotating 
grains should be susceptible to disalignment only during short intervals of slow rotation that 
they called "crossovers"!. Therefore for sufficiently infrequent crossovers suprathermally 
rotating grains will be well aligned with the degree of alignment determined by the time 

1 Radiative torques suggested in Draine & Weingartner (1996) as a means of suprathermal 
rotation are effective only for grains larger than 5 x 10 -6 cm. 

2 Crossovers are due to various grain surface processes that change the direction (in body- 
coordinates) of the systematic torques. 
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between crossovers, the degree of correlation of the direction of grain angular momentum 
before and after a crossover (SM79), and environmental conditions (e.g., magnetic field 
strength B). 

The original calculations of SM79 obtained only marginal correlation of angular 
momentum before and after a crossover, but their analysis disregarded thermal fluctuations 
within the grain material. Lazarian & Draine (1997, henceforth LD97), showed that thermal 
fluctuations are very important, and result in a high degree of correlation for grains larger 
than a critical radius a c , the radius for which the time for internal dissipation of rotational 
kinetic energy is equal to the duration of a crossover. Assuming that the grain relaxation 
is dominated by Barnett relaxation (P79), LD97 found a c 1.5 x 10~ 5 cm, in accord with 
observations that indicated that the dividing line between grains that contribute and those 
that do not contribute to starlight polarization has approximately this value (Kim & Martin 
1995). 

Here we report the discovery that a new effect of thermal fluctuations - which we 
term "thermal flipping" - should lead to alignment of even the small grains with a ^ a c , 
if they rotate suprathermally. However, small grains are observed to not be aligned. We 
argue that this is due to a second effect - "thermal trapping" - which causes small grains to 
rotate thermally a significant fraction of the time, despite systematic torques which would 
otherwise drive suprathermal rotation. 

In §2 we review the role of Barnett fluctuations in the crossover process. In §3 
we discuss how crossovers influence grain alignment, and in §4 we argue that "thermal 
trapping" can account for the observed lack of alignment of small grains. 
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2. Crossovers and Thermal Flipping 



SM79 revealed two basic facts about grain crossovers. First, they showed that in the 
absence of external random torques the direction of grain angular momentum J remains 
constant during a crossover, while the grain itself flips over. Second, they demonstrated that 
in the presence of random torques (arising, for instance, from gaseous bombardment) the 
degree of correlation of the direction of J before and after a crossover is determined by J m i n , 
the minimum value of \J\ during the crossover. As a grain tends to rotate about its axis of 
maximal moment of inertia (henceforth referred to as "the axis of major inertia"), SM79 
showed that the systematic torque components perpendicular to this axis are averaged out 
and the crossover is due to changes in Jy due to the component of the torque parallel to the 
axis.0 Midway through the flipover Ju = and J min = J±. 

The time scale for Barnett relaxation is much shorter than the time between crossovers. 
For finite grain temperatures thermal fluctuations deviate the axis of major inertia from J 
(Lazarian 1994). These deviations are given by the Boltzmann distribution (Lazarian & 
Roberge 1997) which, for an oblate grain (e.g., an axbxb "brick" with b > a) is 



is the kinetic energy, and (3 the angle between the axis of major inertia and J. We define 



where the approximation assumes I z ~ 1.5Ij_, as for an axbxb brick with b/a = y/3. The 

3 Indices || and _L denote components parallel and perpendicular to the axis of major 
inertia. 



f(/3)dP = const • sin/3exp [-E k (P)/kT d ] dp 





(2) 




(3) 
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Barnett relaxation time is (P79) 



t B = 8 x 10 7 a 7 _ 5 (J d /J) 2 sec 



(4) 



where a„ 5 = a/10~ 5 cm. For J > J^, the most probable value of (3 for distribution (JJ) has 
J± = Jsin/3 = Jd, while for J < the most probable value of (3 is ir/2. It follows from 
(H) that during suprathermal rotation (J 2 ^> Jj) the fluctuating component of angular 
momentum perpendicular to the axis of major inertia (J]_) ~ Jj. 



of J||. If t c -C t_B, the Barnett fluctuations can be disregarded during a crossover, and 
J± = const m Jd- The corresponding trajectory is represented by a dashed line in Fig. 1. 
Initially the grain rotates suprathermally with (3 « 0; (3 crosses through 7r/2 during the 
crossover and (3 — >• n as the grain spins up again to suprathermal velocities. 

The condition t c = t B was used in LD97 to obtain a critical grain size a c ~ 1.5 x 10~ 5 cm. 
It was shown that t c < £b for a > a c , and paramagnetic dissipation can achieve an alignment 
of ~ 80% for typical values of the interstellar magnetic field. If paramagnetic aligment were 
suppressed for a < a c this would explain the observed dichotomy in grain alignment: large 
grains are aligned, while small are not. 

What spoils this nice picture is that sufficiently strong thermal fluctuations can enable 
crossovers: fluctuations in (3 span [0, tt] and therefore have a finite probability to flip a grain 
over for an arbitrary value of J. The probability of such fluctuations is small for J 2 3> J^, 
but becomes substantial when | J\ approaches Jd- Indeed, it is obvious from (0) that in the 
latter case the probability of f3 ~ 7r/2 becomes appreciable. LD98 show that the probability 
per unit time of a flipover due to fluctuations is 



Whether the grain trajectory is approximately a straight line in Fig. 1, (a "regular 



SM79 defined the crossover time as t c = 2J±/\j»\ where Jii is the time derivative 




(5) 
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crossover"), or two lines connected by an arc (a "thermal flip") depends on the efficacy of 
the Barnett relaxation. Roughly speaking, thermal flipping happens when t t f £ J / J. If 
J J d the ratio of the flipping and crossover time ttf Ac ~ ^bAc- The latter ratio was found 
in LD97 to be equal to (a/a c ) 13//2 . Therefore flipping was correctly disregarded in LD97 
where only grains larger than a c were considered, but should be accounted for if a < a c . 

The last issue is the problem of multiple flips: a grain with (3 > tt/2 can flip back. 
Thermal flips do not change J. Therefore after a flip (from quadrant (3 < n/2 to quadrant 
(3 > 7i / 2 in Fig. 1) the grain has the same J as before the flip. For J > J d , the thermal 
distribution (M has two most probable values of (3: /3_ = sin _1 (J rf /J), and (3 + = n — (3 + . 
For both (3- and (3 + we have J± = Jd- If we idealize the grain dynamics as consisting 
of systematic torques changing J\\ with J± = const, plus the possibility of instantaneous 
"flips" (at constant J) between f3_ and (3 +) then we can estimate the probability of one 
or more "flips" taking place during a crossover. Let (j)((3)df3 = dt be the probability of 
a flip from (3 to ix — (3 while traversing d/3. The probability of zero flips between and 
(3 is fooift) — ex P[~ Jo 4>(x')dx']. The probability of a "regular crossover" (zero flips) is 
Poo = /oo(tt) = e~ 2a , where 

Similarly, df\o = f'Q <f)d/3 is the probability of one forward flip in the interval d(3, with 
no prior or subsequent flips, and the probability of exactly one forward and zero backward 
flips during the crossover is pio = /io(vr/2) = (1 — e~ 2a )/2. Therefore the probability of one 
or more backward flips is 1 — p 00 — p w = (1 — e~ 2a )/2 — > 1/2 for a a cr . 
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3. Efficacy of Paramagnetic Alignment 

SM79 showed that disalignment of suprathermally rotating grains occurs only during 
crossovers whereas thermally rotating grains undergo randomization all the time. Consider 
a grain subject to random (nonsystematic) torques which provide an excitation rate 
AJ 2 /At, and damping torque —J/td, where td is the rotational damping time. Thermally 
rotating grains have (J 2 ) = (l/2)t rf (AJ 2 /At) = J^ h . This definition of thermally rotating 
grains encompasses grains whose rotation is excited by random H 2 formation, cosmic ray 
bombardment etc. - so long as the associated torques have no systematic component. For 
suprathermally rotating grains J 2 3> Jfh- 

In what follows we roughly estimate the efficacy of grain alignment for t c ^> £b, i-e., 
a < a c . Following P79 we assume that H 2 torques are the dominant spin-up mechanism. 

A crossover requires N ~ J m[n / (AJ Z ) impulse events, where the mean impulse per 
recombination event (see SM79) (AJ Z ) « (2m H a 2 i?/3z/) 1 ^ 2 where v is the number of active 
recombination sites, E is the kinetic energy per newly-formed H 2 , and J min is the minimum 
J reached during the crossover. If iV is multiplied by the sum of mean squared random 
angular momentum impulses ((AJ|) + (AJ 2 )) it gives the mean squared change of J during 
a crossover. Therefore the mean squared change of angle during a flipping-assisted crossover 
is 

^ iV((AJ z 2 ) + (AJ 2 )) ^ ((AJf) + (AJ 2 )) 
{ I J2~ J ■ (AJ) { } 

which differs only by a factor of order unity from the expression for disorientation parameter 
F in SM79, provided that J min is used instead of Jj_. The latter is the major difference 
between the regular crossovers that were described by SM79 and LD97 and our present 
study. SM79 and LD97 dealt with the case for which flipping is negligible and the 
disorientation was mostly happening when J\\ — > due to the action of regular torques, 
in which case J min ~ Jd- As flipping becomes important, J min > Jd is obtained from the 
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condition t tf (J min ) « J min /^- 

Grain alignment is measured by a = (3/2)((cos 2 0) — 1/3) where # is the angle between 
the magnetic field direction and J. Generalizing LD97, 



a A 



3 / arctanye^ff— T 

i + — — , ; 1 

<W V V e cff - 1 , 



1 - A)(7 th 



Jeff 



2Ct b /t r 2.6Ct d /t r ( t 



[l-exp(-(F))] " [l-exp(-(F))] V 1 + t J (9) 
where A = C = 1 in LD97 theory, cr t h is the alignment parameter for thermally rotating 

grains (Roberge & Lazarian 1998), t r is the paramagnetic damping time (Davis & Greenstein 

1951), ib ~ 1.3 + t^) is the mean time back to the last crossover (P79), and ti is the 

resurfacing time. For typical ISM conditions (% = 20 cm~ 3 , B = 5/xG) td/t r ~ 0.05/a_ 5 . 

Expressions (§) and @ (with A = C = 1) were obtained in LD97 assuming that grains 
spend nearly all their time rotating suprathermally, except for brief crossover events with a 
characteristic disorientation parameter (F). We now argue that a significant fraction of the 
small grains do not rotate suprathermally, and an appreciable fraction of crossovers have 
(F) > 1. 



4. Thermal trapping of small grains 

P79 theory of suprathermal rotation did not take into account the "thermal flipping" 
process discussed here. We now argue that thermal flipping will suppress the suprathermal 
rotation of very small grains. 

With the Barnett relaxation time £3 from (|j), the ratio t&/td ~ 
1 x 10~ 4 (15/Td)( Jd/ J) 2 a^_ 5 , where the drag time td is evaluated for a diffuse cloud 
with rin = 30 cm -3 and T = 100K. Thus the timescale for a thermal flip is 

t ti /t d « 10-V 5 (J (i /J) 2 exp{(l/2) [(J/J d ) 2 - l]} (10) 
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showing that thermal flipping is strongly favored for small grains. The critical question now 
is: Can the systematic torques drive small grains to large enough J to suppress the thermal 
flipping, or is the thermal flipping sufficiently rapid to suppress the superthermality? 

Consider a grain with a systematic torque (G — l) 1//2 JthA<2 along the major axis (fixed 
in grain coordinates). The condition J min = J ■ t t f(J min ) becomes 

tti/t d = (J min /J th )/(G-l) 1/2 . (11) 

Thermal flipping causes the systematic torque to randomly change sign in inertial 
coordinates, so that 

(J 2 ) = Jl + (G-l)Jlt tf /(t t{ + t d ) , (12) 
giving a condition for t t f in terms of (J 2 ): 

t*/U = f((JMh) 2 > - ll / \G - ((J/J th ) 2 )} • (13) 



For given a, G, and (Jth/Jd) 2 , ( jl^D an d (0) have either one or three solutions for (J 2 ). 
If (J 2 ) 1 / 2 has multiple solutions J\ < J 2 < J3, the intermediate solution J 2 is unstable: if 
Ji < J < J 2 , then t t f from (|T0|) is smaller than the value required by ([13|), so J — > Ji, 
if J 2 < J < J3, then J — > J 3 . In the former case thermal flipping leads to suppression 
of suprathermal rotation: if the grain enters the region J < J2, then it is trapped with 
J ~ J\ until a fluctuation brings it to J > J2. The timescale for such a fluctuation is 
^trap ~ ^exp[(J 2 /J t h) 2 ]- We refer to this phenomenon as "thermal trapping". 

As an example, consider (Jth/ ' J&) 2 — 5, G — 10 3 and a_ 5 = 0.5. Thermal flipping takes 
place during a crossover at J^ in ~ 5.9J 2 h . If J 2 drops below Jf = 4.5J 2 h , the grain will 
be thermally trapped. For this case thermal flipping will tend to maintain the grain at 
(J 2 ) pa Jf pa 1.02(J t h) 2 , unless thermal fluctuations succeed in getting the grain to J > J 2 , in 
which case thermal flipping is unable to prevent the grain spinning up to a superthermality 
(J/ Jth) 2 ~ G = 10 3 . For this example the thermal trapping time t trap ~ 50^^. During this 
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time paramagnetic alignment will be minimal. Grains that escape the thermal trap become 
suprathermal and align on the timescale for paramagnetic alignment. 

Let rj be the probability per crossover of becoming "thermally trapped" . The fraction 
of grains which are not thermally trapped at any time is x = + ^trap)- 

We can estimate the grain alignment by using (§) and (§) with A = x and 
C = [1 - exp(-(F))]{77 + (1 - V )[l - exp(-(F))]}- 1 . 

During a crossover, the first thermal flip takes place at J ~ J m i n , only a bit larger 
than J2, the thermal trapping boundary. We have seen above that for a < a c , ~ 50% of 
crossovers involve one or more "backward" flips. We do not know what fraction of the 
crossovers end up "thermally trapped" , but we speculate that it could be appreciable, say 

T] ~ 0.1. 

The time between crossovers is of the order of the damping time td (see P79). Returning 
to our example of a grain with a_5 = 0.5, for which we estimated t tTSuP rs 50td, we see 
that the fraction of grains which are not trapped x = 1/(1 + 50rjtd/ib) could be small if 
i] £ 0.1. More detailed studies of the dynamics (Lazarian & Draine 1999) will be required 
to estimate rj, and to provide more reliable estimates of t tvauP , before we can quantitatively 
estimate the degree to which thermal trapping will suppress the alignment of small grains. 

Inspection of Fig. 2 shows that thermal trapping solutions are only found if G is not 
too large (e.g. for G = 10 5 we have no thermal trapping solution in Fig. 2 for a_5 = 0.5). 

Such degrees of suprathermality would follow from variations of accomodation coefficient 
and photoelectric yield, but higher values were obtained in the literature for the case of H2 
torques (P79, LD97). For example, LD97 estimate G = 2 x 10 7 a^ 5 (7/0.2) 2 (10 n cm" 2 /a), 
where a is the surface density of active recombination sites. The values of G ^ 10 4 
required for thermal trapping to be possible for a_5 ~ 0.5 grains would appear to require 
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('j/0.2) 2 /a £ lCT 14 cm 2 . If essentially every surface site is an active chemisorption site, we 
could have a « 10 15 cm~ 2 ; alternatively, it is conceivable that 7 1 for very small grains. 
The latter idea was advocated by Lazarian (1995), who found that oxygen poisoning of 
catalytic sites is exponentially enhanced for grains with a < 10~ 5 cm. Recent experimental 
work (Pironello et al. 1997a,b) suggests that 7 may be much smaller than is usually 
assumed. Moreover, recent research (Lazarian & Efroimsky 1998, Lazarian & Draine 1999) 
shows that faster processes of internal relaxation are possible. These processes should 
enable thermal trapping for larger values of G. 

5. Conclusions 

We have found that "thermal flipping" is a critical element of the dynamics of small 
(a £ 10 -5 cm) grains. If small grains rotate suprathermally, then thermal flipping would 
promote their alignment by suppressing disalignment during "flipping-assisted" crossovers. 
Since small grains are observed to not be aligned, it follows that most must not rotate 
suprathermally. 

One way for small grains to not rotate suprathermally would be for the systematic 
torques from H 2 formation and photoelectric emission to be much smaller than current 
estimates. However, we also find that thermal flipping can result in "thermal trapping", 
whereby rapid thermal flipping can prevent systematic torques from driving small grains to 
suprathermal rotation rates. As a result, at any given time an appreciable fraction of small 
grains are thermally trapped and being disaligned by random processes. 

The thermal trapping effect is of increasing importance for smaller grains, and may 
explain the observed minimal alignment of a £ 5 x 10~ 6 cm dust grains (Kim & Martin 
1995). 
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Fig. 1. — Grain trajectory on the J± - J\\ plane, where J± and J\\ are components of J 
perpendicular or parallel to the grain's principal axis of largest moment of inertia. The solid 
trajectory shows a "thermal flip", while the broken line shows the "regular" crossover which 
would occur in the absence of a thermal flip. 
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Fig. 2.— Upper panel: Eq. flT^) [broken lines, labelled by a_s] and eq. (pTT| ) [solid lines, 
labelled by G]. Dot shows J m i n for flipping-assisted crossover of a_ 5 = 0.5 grain with G = 10 3 . 
Lower panel: Eq. ( [TUP [broken lines, labelled by a_s], and eq. (^3[) [solid lines, labelled by 
G]. Dots show "equilibrium" solutions J\ and J2 for a_5 = 0.5 and G = 10 3 (see text). 



